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ABSTRACT 


Particle  simulations  of  the  modified  two-stream  instability 


demonstrate  strong  electron  acceleration  rather  than  bulk  heating  when 
the  relative  drift  speed  v^  is  below  a  critical  speed  Vc.  A  very 
interesting  nonlinear  autoreaonance  acceleration  process  is  observed 
which  accelerates  the  electrons  much  above  the  phase  speed  of  the 


linearly  unstable  modes.  Simple  criteria  are  presented  that  predict  the 
value  of  Vc  and  the  number  density  of  the  accelerated  electrons. 
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The  modified  two-stream  instability  has  been  studied  rather 
1-3 

extensively  as  an  Important  process  that  transfers  .'on  cross-field 
drift  energy  to  parallel  electron  energy  when  the  electron  ion  drift 
speed  Vj  is  lower  than  the  local  Alfven  speed  v^.  It  is  usually  assumed 
in  the  literature  that  the  electron  energy  transfer  results  in  bulk 
electron  heating*  However,  in  this  letter  we  demonstrate  that  In  the 
low  drift  regime  (i.e.  vd  £  3  Vj,  where  Vj  -  (2T1/m1)1^2  is  ion  thermal 
speed),  the  modified  two-stream  instability  results  in  electron 
acceleration  instead  of  electron  heating  (i.e.  formation  of  suprathermal 
electron  tails).  This  result  has  very  Important  implications  to  many 
physical  situations  such  as  shock  electron  acceleration  *  ,  critical 
ionization  phenomena^  and  the  more  recently  observed  anomalous  glow  of 
the  space  shuttle^.  The  emphasis  in  this  letter  is  on  the  physics  of 
the  tall  formation  as  revealed  by  particle  simulations  rather  than  the 
applications  which  will  be  discussed  elsewhere. 

The  simulations  were  performed  in  one  dimension  in  a  fashion 
similar  to  that  of  McBride  et  al Namely,  the  magnetic  field  (Bq)  was 
tilted  by  an  angle  9  from  the  simulation  axis  (x-axis)  such  that  the 
maximum  linear  growth  of  the  instability  occurs  in  this  direction.  A 
Darwin  electromagnetic  code  with  periodic  boundary  conditions  and  a 
large  system  (128  or  256  cells)  was  used  for  the  computations.  Finite 
size  particles  (50  electrons  and  ions  in  each  cell)  were  assigned  three 
velocity  components  randomly  using  a  quiet  start  method.  The  electrons 
were  treated  as  magnetized  and  the  ions  as  unmagnetized.  The  ions  had 
an  initial  cross-field  drift  relative  to  the  electrons,  ranging  from 
2  to  5  times  v^,  while  the  electrons  were  initially  at  rest.  The  cell 
size  A  was  chosen  as  A  ■  (*e  ■  18  tl1*  ^^ye  length 
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where  vg  -  (2T  /me)*^  is  electron  thermal  speed)  and  the  particle  size 

ax  as  a^/L  “  ^»008  where  L  ■  128A  or  256A  Is  the  system  size.  Other 

initial  parameters  were:  electron  cyclotron  frequency  w,e  »  electron 

plasma  frequency  w  ;  ratio  of  electron  temperature  to  ion  temperature 
pe 

Te/TA  *  2/3;  ion  beta  8  ■  8nnTi/B0^  *  0.01;  mass  ratio  mj/nig  »  400;  the 
angle  9  -  88°  for  the  v^/v^  ■  2  and  9  *  85°  for  v^/vj  »  5  case,  as  given 
by  the  most  unstable  mode. 

The  time  development  of  Fourier  amplitudes  of  the  electric  field 

1/2 

E  (t)  normalized  by  (SitnT  )  is  shown  in  Fig.  1.  (Results  of  the 
*  1 

v^/v^  »  2  run  are  described  unless  otherwise  specified.)  For  this  run 

(L  ■  256A),  approximately  30  modes  are  actually  unstable;  however,  only 

the  most  unstable  modes  m  «  11,  12,  13  are  plotted.  Between  u>LHt  "  5 

and  15,  the  waves  grow  exponentially  and  saturate  at  about  w^t  ~  25. 

At  <»>LHt  ~  20,  the  most  dominant  modes  are  m  *  10  ~  15.  As  will  be  noted 

later,  the  wave  characteristics  change  drastically  between  w^t  •  20  and 

30  due  to  the  wave-particle  Interactions.  At  ~  30,  after 

saturation,,  the  dominant  modes  are  m  ■  6  ~  16  and  m  *  21  ~  26,  while 

after  w^t  ~  35  the  wave  amplitudes  decrease  slowly.  To  show 

contributions  of  the  electric  fields  other  than  m  -  11  ~  13  at  later 

2 

times,  electric  field  energy  e  *  <E  >/ 8*  is  also  plotted  in  Fig.  1. 

h  x 

The  linear  theory  of  the  modified  two-stream  instability  is 

3 

described  in  the  electrostatic  limit  by  the  dispersion  equation 

DE8  5  1  +  <2“p./k“>  I*  +  <0Z«0> 

♦  <2»21/k2v2)  (1  +  -  0, 


(1) 


4 


where  <;Q  -  w/k(ve,  Cj  ■  (w  -  kv^sln  9)/kvt>  k(  «  k  cos  9,  Z(c)  Is  the 

plasma  dispersion  function  and  lg(A)  is  the  zero-th  order  modified 

2  2  2  2 

Bessel  function  with  K  *  k  sin  9/2^^.  Hie  measured  frequency  and 

growth  rate  in  our  simulations,  scaled  by  the  lower-hybrid  frequency 

“LH»  are  wr  ~  °*3  “lh*  y  ~  °‘1  WLH  ^or  t*'e  most  unstable  mode,  and  the 

1/2 

corresponding  wavenumber  k  satisfies  kv. /w. „  -  kv  /u>  (T  /T  >  ~  0.4  , 

i  LH  e  ce  i  e-  ’ 

all  of  which  agree  well  with  the  linear  theory. 

The  distribution  function  of  electron  velocities  parallel  to  the 
magnetic  field  at  the  end  of  the  run  is  shown  in  Fig.  2.  Superimposed 
on  them  are  the  corresponding  ion  perpendicular  velocity  distributions 
f(v^)  ,  where  v  =  v^/cos  9  .  A  well-developed  high  energy  electron 
tail  is  formed  with  the  ma.  imum  parallel  speed  of  vjmax  ~  ^ve  while 
little  change  occurs  in  the  main  part  of  electrons,  which  retains  its 
Maxwellian  shape  centered  at  v#  »  0  .  On  the  other  hand,  the  initial 
drifting  Maxwellian  distribution  of  ions  centered  at 
v^  ■  4v^  evolves  into  double-peaked  distribution  which  suggests  strong 
ion  trapping. 

The  p<**-.r=  space  distributions  of  the  electrons  and  ions  are  shown 
in  Fig.  3.  Plotted  here  are  the  electron  parallel  velocities  (left 
column)  and  ion  perpendicular  velocities  (right  column)  versus  their 
position  along  the  x  direction.  Distributions  in  other  directions  do 
not  show  significant  changes  except  for  a  slight  electron  ExB  drift  and 
are  not  displayed  here.  Changes  in  the  distributions  begin  to  occur  at 
~  15  which  is  the  start  of  the  nonlinear  stage.  Between  w^t  ■  20 
and  30,  electrons  are  rapidly  accelerated  along  the  magnetic  field  while 
iona  become  trapped  by  the  waves.  The  electron  phase  space  at  w^t  " 
(Fig.  3b)  shows  13  or  14  tails,  about  the  same  number  as  the  mode  number 
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(12)  of  the  linearly  most  unstable  wave.  After  u>^Ht  ~  30,  the  electrons 

In  the  tails  are  slowly  and  raonotonically  accelerated  without  a 

significant  bouncing  motion,  and  they  finally  become  untrapped  and 

detached  fiom  the  main  body  of  electrons  to  form  a  high  energy  electron 

beam  parallel  to  the  magnetic  field  (see,  Fig.  3c-e).  The  fastest 

2 

electrons  have  energies  mv  12  -  50  T  which  Is  100  times  the  Initial 

#  eo 

electron  energy.  Since  the  population  of  the  electron  beam  Is 
approximately  8%  of  the  total  electrons  and  their  average  energy  is  36 

Teo  or  six  times  as  much  as  the  Initial  Ion  drift  energy,  the  electron 

beam  carries  60%  of  the  total  electron  energy  parallel  to  the  magnetic 

field  at  the  end  of  the  run.  The  saturation  mechanism  of  the 

Instability  Is  ion  trapping,  as  seen  in  Fig.  3.  The  saturation 

2  2 

amplitude  Is.  calculated  from  IF.^I  ~  k  ~  (m^/2)! (v^-w/k)  +Vj/2]. 

The  theoretical  saturation  amplitude  due  to  ion  trapping  is  thus 

(E//8iinTt)th  ~  0.14  and  the  amplitude  measured  in  the  simulation  is 

(E//8nnT  )  .  ~  0.11  for  the  v./v.  *  2  run. 

i  obs  ai 

The  evolution  of  the  wave  characteristics  and  the  electron  tall 
formation  are  strongly  coupled.  Shown  in  Fig.  4  Is  the  change  in  the 
parallel  wave  phase  speed  v  ^  *  w/k^  f°r  Fourier  modes  m  ■  10  ~  13  and 
24.  (The  phase  speed  was  calculated  from  the  measured  frequency  u>.] 

The  two  solid  lines  represent  the  maximum  speed  of  the  electrons 
parallel  to  the  magnetic  field  (denoted  el.  max)  and  the  speed  of  the 
fastest  1%  electrons  on  the  tail  (denoted  el.  1%).  As  the  linear  theory 
predicts,  the  phase  speeds  of  the  unstable  modes  are  all  less  than  1.5 
v€  within  the  slope  of  the  electron  distribution  before  the  instability 
seta  in.  This  means  that  the  electrons  can  resonantly  Interact  with  the 
waves.  During  the  nonlinear  growth  of  the  instability  (ojyjt  ■  20  ~  30), 
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the  phase  speed  of  the  wave  always  stays  at  the  t !.•  •>'  du-  rapidly 

rising  electron  tails  of  Fig.  3.  Ourtng  the  same  Interval  the  phase 

speed  increases  (l.e.,  there  is  nonlinear  frequency  shtft)  closely 

followed  by  an  Increase  In  the  velocities  of  the  electrons  in  the 

tall.  When  the  phase  speed  reaches  the  drift  speed,  v  /cos  9  ■  4v 

d  e 

(projected  on  the  parallel  direction),  it  stops  increasing  and  stays 
close  to  the  drift  speed  for  w^t  "  30  ~  50.  After  u^t  ~  40,  slowly 
growing  waves  with  the  mode  numbers  around  8  and  24  dominate  and 
continue  to  slowly  accelerate  the  electrons.  By  increasing  the  ratio 
vd/vi  to  larger  than  3,  we  observed  a  gradual  transition  from  tail 
formation  to  bulk  heating.  Our  results  were  similar  to  Ref.  3,  for  the 
case  Vj/Vj  “  5,  in  which  bulk  electron  heating  occurs. 

The  observed  critical  speed  V£  ■*  3v^  below  which  electron 
acceleration  occurs.  Is  consistent  with  considerations  concerning 
whether  electrons  or  ions  are  trapped  first  *  .  The  ratio  of  the 
electron  (*e)  to  Ion  (*j)  trapping  potential  is  given  by 

V*e  '  <V*.><Vk>2  <'d  -  V2/’ph2  (2> 

where  v^  is  the  phase  speed  of  the  most  unstable  wave.  For  ®j/*e  >  1 
electron  trapping  occurs  first  and  saturates  the  instability  resulting 
in  bulk  electron  heating.  In  the  opposite  case  <  1),  saturation 

occurs  by  ion  trapping  and  only  a  fraction  of  the  electrons  gain 
energy.  From  (l)  and  the  inequality  <  l,  we  find  the  condition 

for  tail  formation  as 

v. 

U.b  <  —  <  3 
vi  ~ 


(3) 
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The  lower  inequality  comes  from  the  linear  threshold.  Our  runs  for 
v^'i  “  3.5  a°d  "’•Q  were  characterized  by  flat  topped  electron 
distributions  and  absence  of  tails,  which  are  consistent  with  the 
inequality  (3). 

The  observed  fraction  of  electrons  being  accelerated  is  almost 
equal  to  the  fraction  of  electrons,  that  initially  had  velocity  larger 
than  the  phase  speed  v  ^  =  “>/k  B  corresponding  to  the  most  unstable 
mode.  For  example,  for  the  case  v^/v^  *  2  the  above  argument  gives  9Z, 
while  the  observed  value  was  8%.  It  is  interesting  to  note  that  in  this 
case  the  electron  trapping  width  was  6vT  »  1.5  vg,  comparable  to  the 
parallel  wave  phase  velocity.  The  lack  of  trapped  orbits  in  the 
electron  phase  space  Indicates  that  the  autocorrelation  time  was  shorter 
than  the  trapping  time.  The  simulation  results  indicate  that  a 
nonlinear  autoresonance  process  operates  accelerating  electrons  to  their 
final  energy.  While  the  trapped  ion  dynamics  control  the  saturation  of 
the  linearly  unstable  modes,  the  tall  electrons  (Fig.  3)  control  the 
nonlinear  frequency  shift,  which  results  in  an  increase  in  the  parallel 
phase  speed  shown  in  Fig.  4.  This  result  is  supported  by  our  frequency 
spectral  analysis  which  shows  no  frequency  shift  for  the  v^/vj  *  5  case, 
contrary  to  the  v^  <  Vc  case. 

In  summary  we  have  demonstrated  that  for  v^  <  V£  *  3v^,  the 
modified  two-stream  instability  results  in  electron  acceleration  rather 
than  the  conventionally  assumed  electron  heating.  In  this  case  ion 
trapping  controls  the  instability  saturation,  while  a  small  fraction  of 
electrons  are  accelerated  in  the  direction  parallel  to  the  magnetic 
field  by  a  systematic  nonlinear  increase  in  the  phase  velocity.  The 
simulations  shown  here  were  1-D.  Preliminary  2-D  simulations,  which 
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allow  for  a  symmetric  growth  of  the  fields  to  the  magnetic  field  show 
development  of  similar  electron  tails.  In  this  case,  however,  the  tails 
are  on  both  sides  of  the  eLeciron  distribution.  These  results  along 
with  a  detailed  analysis  of  the  autoacceleration  process  is  under  study 
and  will  be  discussed  elsewhere. 
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Figure  Captions 


Fig.  1  Development  of  three  most  unstable  Fourier  amplitudes  (labelled 

by  node  number)  of  the  electric  field  normalized  by 

1/2 

(8unT^)  where  n  and  are  initial  density  and  ion 

temperature.  Also  plotted  is  electric  field  energy 

t  =  <E2>/8it. 

E  x 

Fig.  2  Electron  parallel  (solid)  and  ion  perpendicular  (dashed) 

velocity  distribution  at  u>  t  =  60.  The  ordinate  is  in  a 

LH 

linear  scale. 

Fig.  3  Phase  space  distributions  of  electrons  (v./v  vs  x)  (left)  and 

if  e 

ions  (v^/v  vs  x)  (right)  at  w^t  =  20,  30,  40,  50  and  60  from 
top  to  bottom. 

Fig.  4  Evolution  of  the  wave  phase  speed  parallel  to  the  magnetic 

field  for  modes  m  *  10  ~  13  and  24.  Top  two  lines  represent 
the  maximum  electron  speed  and  the  speed  of  the  fastest  (1%) 
electrons  respectively. 


